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Abstract 

We study static fermion bags in the 1 + 1 dimensional Gross-Neveu and Nambu- 
Jona-Lasinio models. It has been known, from the work of Dashen, Hasslacher and Neveu 
(DHN), followed by Shei's work, in the 1970's, that the self-consistent static fermion bags in 
these models are reflectionless. The works of DHN and of Shei were based on inverse scattering 
theory. Several years ago, we offered an alternative argument to establish the reflectionless 
nature of these fermion bags, which was based on analysis of the spatial asymptotic behavior 
of the resolvent of the Dirac operator in the background of a static bag, subjected to the 
appropriate boundary conditions. We also calculated the masses of fermion bags based on the 
resolvent and the Gelfand-Dikii identity. Based on arguments taken from a certain generalized 
one dimensional supersymmetric quantum mechanics, which underlies the spectral theory of 
these Dirac operators, we now realize that our analysis of the asymptotic behavior of the resolvent 
was incomplete. We offer here a critique of our asymptotic argument. 
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1 Introduction 



Many years ago, Dashen, Hasslacher and Neveu (DHN) [T[], and following them Shei 
0, used inverse scattering analysis to find static fermion-bag ^ soliton solutions 
to the large-iV saddle point equations of the Gross-Neveu (GN) [0 and of the 1 + 1 
dimensional, multi-flavor Nambu-Jona-Lasinio (NJL) |^ models. In the GN model, 
with its discrete chiral symmetry, a topological soliton, the so called Callan-Coleman- 
Gross-Zee (CCGZ) kink 0, was discovered prior to the work of DHN. 

One version of writing the action of the 1 + 1 dimensional NJL model is 

S = Jd^xlY^^a [z^-((7 + Z7r75)] ^Pa ' ^ {cT^ + 7^^)] , (l-l) 



.a=l 



where the ipa{a- = 1, • • • , N) are N flavors of massless Dirac fermions, with Yukawa 
couplings to the scalar and pseudoscalar auxiliary flelds cr(x), 7r(x)[|. 

The remarkable discovery DHN made was that all these static bag conflgurations 
were reflectionless. More precisely, the static cr(x) and tt{x) conflgurations, that solve 
the saddle point equations of the NJL model, are such that the the Dirac equation 

i0 — — i7r(a;)75 'i}j{x) = (1.2) 

in these backgrounds has scattering solutions, whose reflection amplitudes at momen- 
tum k vanish identically for all values of k. In other words, a fermion wave packet 
impinging on one side of the potential well cr(x) +i7r(x)75, will be totally transmitted 
through the well (up to phase shifts, of course). 

We note in passing that besides their role in soliton theory reflectionless 
potentials appear in other diverse areas of theoretical physics |T^, [IT], [T2|. For a 
review, which discusses reflectionless potentials (among other things) in the context 
of supersymmetric quantum mechanics, see ||T3| . 

Since the works of DHN and of Shei, these fermion bags were discussed in the 
literature several other times, using alternative methods [Q. For a recent review on 



^The fermion bag solitons in these models arise, as is well known, at th e le vel of the effective 



action, after integrating the fermions out, and not at the level of the action (1.1 ) 
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these and related matters, see |T5[. Very recently, static chiral fermion bag solitons 



IGf in a 1 + 1 dimensional model, as well as non-cliiral (real scalar) fermion bag 



solitons |T^, were discussed, in which the scalar field that couples to the fermions 
was dynamical already at the classical level (unlike the auxiliary fields a and vr in 

O))- 

In many of these treatments, one solves the variational, saddle point equations 
by performing mode summations over energies and phase shifts. An alternative to 
such summations is to solve the saddle point equations by manipulating the resolvent 
of the Dirac operator as a whole, with the help of simple tools from Sturm-Liouville 
operator theory. The resolvent of the Dirac operator takes care of mode summation 
automatically. 

Some time ago, one of us had developed such an alternative to the inverse scat- 
tering method, which was based on the Gel'fand-Dikii (GD) identity ||T^ (an identity 
obeyed by the diagonal resolvent of one-dimensional Schrodinger operators)^, to study 



fermion bags in the GN model [ 19 1 as well as other problems [BO]. That method was 



later applied by us to study fermion bags in the NJL model and in the massive 
GN model [^. Similar ideas were later used in |^ to calculate the free energy of 
inhomogeneous superconductors. 

Application of this method in |]l9l and in [^] reproduced the static bag results of 
DHN and of Shei in what seems to be a simpler manner than in the inverse scattering 



formalism. In we followed the method introduced in ||19| , |2(| , and simply wrote 
down an efficient, parameter dependent, ansatz for the diagonal resolvent of the Dirac 
operator in a static ct{x), t{{x) background. Construction of that ansatz was based 
on simple dimensional analysis, and on the Gelfand-Dikii identity. Nowhere in the 
construction, did we use the theory of refiectionless potentials. With the help of that 
ansatz, we were able to reproduce in Shei's inverse scattering results, in a similar 
manner to the reproduction of DHN's results in p!9| . 

In addition to rederivation of bag profiles, masses and quantum numbers found 



-For a simple derivation of the GD identity, see l2(j, 21 
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by DHN and Shei, we tried in |^ to explain the reflectionless property of the static 
background in simple terms, by studying the expectation value of the fermion current 
{j^{x)) = {i/j{x)'y^'ilj{x)) in a given background cr(x) + i'y^n{x), at spatial asymptotics. 
However, after careful reexamination, we now realize that the part of the analysis 
pl| on the reflectionless nature of the background was incomplete.'^ We realized 



m 



this with the help of a certain version of generalized one dimensional supersymmetric 



(SUSY) quantum mechanics p4|, that underlies the spectral theory of the Dirac 
operator in (|1.2|). 



This paper offers critique of our asymptotic argument from The rest of the 
results in [^, namely, bag profiles etc., remain intact, and will not be discussed here. 

Before discussing this issue in detail, and in order to set our notations, let us recall 
some basic facts about dynamics of the NJL model: 

The partition function associated with ( |1 . 1|) isQ 

Z = J VaVTrVipVip expi J rf^x [i^ - (a + ^75)] - ^ (o-^ + vr^) | (1.3) 

Integrating over the grassmannian variables leads to Z = J VaVn exp{iS'e//[cr, vr]} 
where the bare effective action is 



S, 



1 



2g' 



d X 



a 



IT 



iN Tr log i0 — (a 



;i-4) 



and the trace is taken over both functional and Dirac indices. 

This theory has been studied in the limit N 00 with Ng'^ held fixed 0]. In this 
limit ( |1.3| ) is governed by saddle points of ( |1.4| ) and the small fluctuations around 
them. The most general saddle point condition reads 



6S^ 



eff 



a {x, t) 



5a (x, t) 



+ iN tr 



1 

(cr + Z7r75) 



\x,t) 



■^Wc thank R. Jaffe and N. Graham for useful correspondence on this point. 

^From this point to the end of this paper flavor indices are usually suppressed. Thus i 



should 



N 



N 



be understood as i '^^^ ^a$ipa- Similarly ^pTij} stands for ipa^ipa, where T = 1,75. 



a=l 
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- (a + Z7r75) 



0. (1.5) 



In particular, the non-perturbative vacuum of ( |1.1| ) is governed by the simplest 
large N saddle points of the path integral associated with it, where the composite 
scalar operator ipip and the pseudoscalar operator iil!'y5ip develop space-time indepen- 
dent expectation values. 

These saddle points are extrema of the effective potential Veff associated with 
(|1.1|), namely, the value of —Sejf for space-time independent a,7i configurations per 
unit time per unit length. The effective potential Ve/f depends only on the combina- 

eff 

of p at p = m 7^ that is fixed by the (bare) gap equation | 

d^k 1 



tion p^ = o"^ + TT^ as a result of chiral symmetry. Veff has a minimum as a function 



_,„ + ,A.,nr/p-^^ = (1,6) 
which yields the dynamical mass 

m = Ae~^^^. (1.7) 

Here A is an ultraviolet cutoff. The mass m must be a renormalization group invariant. 
Thus, the model is asymptotically free. We can get rid of the cutoff at the price of in- 
troducing an arbitrary renormalization scale p. The renormalized coupling g^^ (p) and 

1 J 

the cut-off dependent bare coupling are then related through A e '^a^i^) = pe 

in a convention where Ng\{m) = K Trading the dimensionless coupling g\ for 

the dynamical mass scale m represents the well known phenomenon of dimensional 

transmutation. 

The vacuum manifold of (|1 . 1|) is therefore a circle p = m in the a, tt plane, and 
the equivalent vacua are parametrized by the chiral angle 9 = arctan^. Therefore, 
small fluctuations of the Dirac fields around the vacuum manifold develop dynamical 
chiral mass m expliO'j^) . 

Note in passing that the massless fluctuations of 6 along the vacuum manifold 



decouple from the spectrum |^5[ so that the axial U{1) symmetry does not break 



dynamically in this two dimensional model in accordance with the Coleman- 
Mermin- Wagner theorem. 

Non-trivial excitations of the vacuum, on the other hand, are described semi- 



classically by large N saddle points of the path integral over ( |1 . 1|) at which a and vr 
develop space-time dependent expectation values]^ ^. These expectation values are 
the space-time dependent solution of ( p..5| ). Saddle points of this type are important 
also in discussing the large order behavior ||28|, |29[| of the expansion of the path 
integral over 

These saddle points describe sectors of ( |1.1|) that include scattering states of the 
(dynamically massive) fermions in (p..!]), as well as a rich collection of bound states 
thereof. 

These bound states result from the strong infrared interactions, which polarize the 
vacuum inhomogeneously, causing the composite scalar ipilj and pseudoscalar iip'^r^ip 
fileds to form finite action space-time dependent condensates. These condensates are 
stable because of the binding energy released by the trapped fermions and therefore 
cannot form without such binding. This description agrees with the general physical 



picture drawn in We may regard these condensates as one dimensional chiral 
bags H that trap the original fermions ( "quarks" ) into stable finite action extended 
entities ("hadrons"). 

If we set 71 (x) in ( |1 . 1|) to be identically zero, we recover the Gross- Neveu model, 
defined by 

Sgn = J d^x i^ij[t0-a]^-^Y (1.8) 

In spite of their similarities, these two field theories are quite different, as is well- 
known from the field theoretic literature of the seventies. The crucial difference is 
that the Gross- Neveu model possesses a discrete symmetry, a —a, rather than 
the continuous axial U{1) symmetry a + 27577 e~'^"'^''^{a + 27577) in the NJL model 



TTT|) . This discrete symmetry is dynamically broken by the non-perturbative vacuum. 



and thus there is a kink solution |^, |I|, |T^, the CCGZ kink mentioned above, (t{x) = 
mtanh(ma;), interpolating between ±m at x = ±00 respectively. Therefore, topology 
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insures the stability of these kinks. 

In contrast, the NJL model, with its continuous symmetry, does not have a topo- 
logically stable soliton solution. The solitons arising in the NJL model can only be 
stabilized by binding fermions, namely, stability of fermion bags in the NJL model is 
not due to topology, but to dynamics. 

The rest of the paper is organized as follows: In Section 2 we review the results 



of pj]. We study the resolvent of the Dirac operator in a given static a{x) +275'7r(x) 
background. The Dirac equation in any such background has special properties. 
In fact, we show that it is equivalent to a pair of two isospectral Sturm-Liouville 
equations in one dimension, which generalize the well known one-dimensional super- 
symmetric quantum mechanics. We use this generalized supersymmetry to express all 
four entries of the space-diagonal Dirac resolvent (i.e., the resolvent evaluated at co- 
incident spatial coordinates) in terms of a single function. As a result, we can prove 
that each frequency mode of the spatial current {il!{x)'~f^tp{x)) vanishes identically. 



contrary to the argument we made in |21|. The findings of Section 2 are then used 



in Section 3 to simplify the saddle point equations ( pTSl) . We then study the spa- 
tial asymptotic behavior of the simplified equations. We use the spatial asymptotic 
expression of the resolvent of the Dirac operator (summarized in the Appendinx) to 
generate an asymptotic expansion of the quantities and ^^i^^, evaluated on a 

static background ((t(x), 7r(x)) (consistent with the physical boundary conditions at 
spatial infinity). We prove that these asymptotic expansions vanish term by term to 
any power in 1/x, for any static a{x) and 7r(x) that are consistent with the physical 
boundary conditions, and not just for refiectionless backgrounds, as we have claimed 
PT| . In the Appendix we recall the asymptotic behavior of resolvents of Sturm- 



m 



Liouville operators and use them to derive the asymptotic behavior of the resolvent 
of the Dirac operator in a static bag background. 
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2 Resolvent of the Dirac Operator With Static 
Background Fields 



As was explained in the introduction, we are interested in static space dependent 
solutions of the extremum condition on Sf-g. To this end we need to invert the Dirac 
operator 

D = i$-{a{x)+m{x)-i5) (2.1) 

in a given background of static field configurations a{x) and 7r(x). In particular, 
we have to find the diagonal resolvent of ( ^.11 ) in that background. We stress that 
inverting ( |2.1| ) has nothing to do with the large approximation, and consequently 
our results in this section are valid for any value of A^. For example, our results may 
be of use in generalizations of supersymmetric quantum mechanics. 

For the usual physical reasons, we set boundary conditions on our static back- 
ground fields such that (t{x) and 7r(x) start from a point on the vacuum manifold 
cr^ + vr^ = at x = — oo, wander around in the a — ir plane, and then relax back 
to another point on the vacuum manifold at x = +oo. Thus, we must have the 
asymptotic behavior 

a — > mcos9± , a' — > 

TT — , msme± , vr' — , (2.2) 

I'— +±oo a;— »itoo 

where 6± are the asymptotic chiral alignment angles. Only the difference 6'+ — 6^ is 
meaningful, of course, and henceforth we use the axial U{1) symmetry to set 6^ = 0, 
such that cr(— oo) = m and 7r(— oo) = 0. We also omit the subscript from 6^ and 
denote it simply by 6 from now on. As typical of solitonic configurations, we expect, 
that (j{x) and 7r(x) tend to their asymptotic boundary values (|2.2|) on the vacuum 
manifold at an exponential rate which is determined, essentially, by the mass gap m 
of the model. It is in the background of such fields that we wish to invert (|2.1| ). 
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In this paper we use the Majorana representation 

7^^ = (72 , 7^ = io"3 and 7^ = —7^7^ = o"i (2.3) 

for 7 matrices. In this representation becomes 

(—dx — o —ioo — m \ I —Q —iuo — ivr \ 
(2.4) 
ioj — in dx — (T / \ iu — iTC — / 

where we introduced the pair of adjoint operators 

Q = a{x) + dx, = a{x) - . (2.5) 



-il5{x-y) (2.6) 



(To obtain ( p. 41) , we have naturally transformed i0— ((t(x) +i7r(x)75) to the uj plane, 
since the background fields a{x),iT{x) are static.) 
Inverting ( p.4|) is achieved by solving 

^ —Q —iuj — m{x) \ I a{x,y) b{x,y) 

\iuj-iTi{x) -Q^ J \ c{x,y) d{x,y) 

for the Green's function of (|2.4|) in a given background a{x),TT{x). By dimensional 
analysis, we see that the quantities a, b, c and d are dimensionless. 

2.1 Generalized "Supersymmetry" in a Chiral Bag Back- 
ground 



Interestingly, the spectral theory of the Dirac operator (|2.4| ) is underlined by a certain 



generalized one dimensional supersymmetric quantum mechanics ||2^. This general- 
ized supersymmetry is very helpful in simplifying various calculations involving the 
Dirac operator and its resolvent. In the remaining part of this section, we review the 
discussion in p4 . 



The diagonal elements a{x,y), d{x,y) in ( |2.6| ) may be expressed in term of the 
off-diagonal elements as 

aix, y) = --Q^c{x, y) , d{x, y) = — ■ — ^Qb{x, y) (2.7) 



which in turn satisfy the second order partial differential equations 
1 



Q 



t. 



uj + vrfx 



-Q - {ijj - 7r(a;)) 



h{x,y) 



Q 



dxb{x, y) 


+ 


UJ + 7r(x) 





cij + t:[x) 



Kx,y) 

oj + 7r(x) 



5{x - y) 



1 



-a 



UJ — Tx{x) 

dxc{x, y) 
OJ — 7r(a;) 



Qt _ + 7r(a;)) 



c{x,y) 



+ 



2 , ajxYix) 
UJ — ^{x) 



c{x,y) 
UJ — 7r(x) 



-5{x - y) 
(2.8; 



Thus, ?/) and — c(x, ?/) are simply the Green's functions of the corresponding 
second order Sturm-Liouville operators^ 

h{x) 



Lb{uj)b{x) = -dx 



LJuj)c(x) = -dx 



in (|2.8|), namely, 



dxb{x) 


+ 


UJ + 7r(x) 


dxc{x) 


+ 


UJ — 7^{x) 



a{xf + 7r(x)2 - a'(x) - a;^ + ^(^KM 



a{xf + 7r(a;)^ + a'(a;) - uj^ + ^(^KM 

a; — TT X 



c<j + 7r(x) 

c(x) 
— 7r(x) 
(2.9) 



b{x,y) 
c{x,y) 



9 (.T - y) h2{x)hi{y) + 9 {y - x) h2{y)hi{x) 

9{x-y) C2{x)ci{y) + 9 {y - x) C2{y)ci{x) 
W, 



(2.10) 



Here {hi{x),h2{x)} and {ci(x), C2(x)} are pairs of independent fundamental solutions 
of the two equations L^b^x) = and Lcc{x) = 0, subjected to the boundary conditions 

6i(x),ci(x) <c(fc)e-^^^ , 



&2(x),C2(x) — . A,,(A;)(2)e^'=-' (2.11) 

x— >+oo 



^Notc that OJ plays here a dual role: in addition to its role as the spectral parameter (the 



terms in ( |2.9| )), it also appears as a parameter in the definition of these operators- hence the explicit 
Lu dependence in our notations for these operators in (|]^). However, in order to avoid notational 
cluttering, from now on we will denote these operators simply as Lb and Lc- 
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with some possibly k dependent coefficients A^f^l{k), A^^l{k) and withQ 



k = Vu^^^, ImA;>0. (2.12) 

The purpose of introducing the (yet unspecified) coefficients A^l{k), A^^l{k) will be- 
come clear following Eqs. ( p.l5|) and (|2.16|) . The boundary conditions (|2.11| ) are 



consistent, of course, with the asymptotic behavior ( p.2| ) of a and tt due to which 
both Lh and tend to a free particle hamiltonian [—d"^ + — u"^] as x ^ ±00. 
The wronskians of these pairs of solutions are 

b2{x)b[{x) - bi{x)b'2{x) 



Wt{k) 
W,{k) 



uj + 7r(x) 

C2{x)c^{x) - Ci(x)c2(x) 

uo — 7r(x) 

(2.13) 



As is well known, Wi,{k) and Wc{k) are independent of x. 

Note in passing that the canonical asymptotic behavior assumed in the scattering 
theory of the operators Li, and corresponds to setting A^^l = A^^l = 1 in (|2.11|) . 
Thus, the wronskians in ( |2.13| ) are not the canonical wronskians used in scattering 
theory. As is well known in the literature [Q], the canon^ca/ wronskians are proportional 
(with a k independent coefficient) to k/t{k), where t{k) is the transmission amplitude 
of the corresponding operator Lf, or Lc- Thus, on top of the well-known features 
of t{k), the wronskians in (|2.13|) will have additional spurious /c-dependence coming 
from the amplitudes A^f^l{k), A^^l{k) in ( |2.11| ). 



Substituting the expressions ( p.lO|) for the off-diagonal entries b{x, y) and c(x, y) 
into ( |2.7D , we obtain the appropriate expressions for the diagonal entries a{x,y) and 
d{x,y). We do not bother to write these expressions here. It is useful however to 
note, that despite the ^^.'s in the Q operators in ( |2.7| ), that act on the step functions 
in ( p.lO| ), neither a{x,y) nor d{x,y) contain pieces proportional to 6{x — y) . Such 



pieces cancel one another due to the symmetry of (|2.10|) under x ^ y. 

^We see that if Im/c > 0, 61 and ci decay exponentially to the left, and 62 and C2 decay to the 
right. Thus, if Im/c > 0, both b{x, y) and c(a;, y) decay as \x — y\ tends to infinity. 
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We will now prove that the spectra of the operators and Lc are essentially the 
same. Our proof is based on the fact that we can factorize the eigenvalue equations 
Lhb{x) = and Lcc{x) = as 



UJ — Tlix) UJ + TT[x] 



■ Q TT c 



uj + 7r(x) uj — 7r(x) 

(2.14) 



as should be clear from (|2^ ) and (|2.9|) . 

The factorized equations ( p.l4| ) suggest the following map between their solutions. 
Indeed, given that Li,b{x) = 0, then clearly 

is a solution of Lcc{x) = 0. Similarly, if Lcc{x) = 0, then 

b{x) = ^-—Q^c{x) (2.16) 

UJ — TC[X) 

solves Lbb{x) = 0. 

Thus, in particular, given a pair {&i(x), 62(2^)} of independent fundamental so- 
lutions of Lbb{x) = 0, we can obtain from it a pair {ci(x), C2(x)} of independent 
fundamental solutions of Lcc{x) = by using (|2.15|) , and vice versa. Therefore, with 
no loss of generality, we henceforth assume, that the two pairs of independent funda- 
mental solutions {61 (x), 62(3;)} and {ci(x), C2(x)}, are related by ( p.l5|) and (|2.16|) . 



The coefficients A^f^l{k), Al^l{k) in ( ^.11| ) are to be adjusted according to ( |2.15| ) 



and (|2.16|) , and this was the purpose of introducing them in the ffist place. 



Thus, with no loss of generality, we may make the standard choice 

= = 1 (2.17) 



in ( p.ll|) . The coefficients A^^^ are then determined by (|2.15| ): 

(1) _ a(-oo) -ifc 

c — — 7 ^ 

TT^ — OO) + U 
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= 44±1^. (2.18) 

" 7r(oo)+u; ^ ^ 

We note that these h{x) ^ c{x) mappings can break only if 

Q6 = or Q^c = 0, (2.19) 



for h{x) or c(x) that solve ( p. 141) . Do such solutions exist? Let us assume, for example, 



that Qh = Q and that L56 = 0. From the first equation in (|2.14| ) (or in (|2.8|) ), we 
see that this is possible if and only if ± 7r(x) = 0, which clearly cannot hold if 
(9^7r(x) 7^ 0. A similar argument holds for Q''" c = 0. Thus, if dx'^{x) 7^ 0, the 
mappings ( p.l5| ) and ( |2.16| ) are one-to-one. In particular, a bound state in implies 
a bound state in Lc (at the same energy) and vice-versa. 

An interesting related result concerns the wronskians Wh and Wc- From ( |2.13| ), 
and from ( p.l5 ) and ( p.l6 ) it follows immediately that for pairs of independent fun- 



damental solutions {&i(a;), 62(3^)} and {ci(x), C2(x)} we have 

= = C162 - C261 = ■ — = Wb . 2.20 

UJ — TC[X) 00 + 7r[x) 

The wronskians of pairs of independent fundamental solutions of and Lc, which 
are related via (|2.15|) and (|2.16| ) are equal! 



To summarize, if dx'n'lx) 7^ 0, and Lc have the same set of energy eigenvalues 
and their eigenfunctions are in one-to-one correspondence. 

If, however, vr =const., then we are back to the familiar "supersymmetric" factor- 
ization 

Qb= (u^ -7T^)b, QQ^ c= {u^ -n^)c, (2.21) 

and mappings 

c{x) = -^—Q b{x) , b{x) = -^—Q^ c{x) . (2.22) 

UJ -\- Tl UJ — IT 

As is well known from the literature on supersymmetric quantum mechanics, the 
mappings ( |2.22|) break down if either Qb = or Q^c = 0, in which case the two 
operators Q^Q and are isospectral, but only up to a "zero-mode" (or rather, an 
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Lo^ = TT^ mode), which belongs to the spectrum of only one of the operator^. The 
case 7r(x) = brings us back to the GN model. Supersymmetric quantum mechanical 
considerations were quite useful in the study of fermion bags in [^. 

The "Witten index" associated with the pair of isospectral operators and L^, 
is always null for backgrounds in which dxTT{x) ^ 0, since they are absolutely isospec- 
tral, and not only up to zero modes. There is no interesting topology associated with 
spectral mismatches of Lb and Lc- This is not surprising at all, since, as we have 
already stressed in the introduction, the NJL model, with its continuous axial sym- 
metry, does not support topological solitons. This is in contrast to the GN model, for 
which TT = 0, which contains topological kinks, whose topological charge is essentially 
the Witten index of the pair of operators ( |2.21| ). 

We note in passing that isospectrality of and which we have just proved, 
is consistent with the 75 symmetry of the system of equations in (|2.6|) , which relates 
the resolvent of D with that of = —'-y^D'-y^. Due to this symmetry, we can map the 
pair of equations Lbb{x, y) = 6{x — y) and Lcc{x, y) = —6{x — y) (Eqs. (p.8|)) on each 
other by 

b{x,y) ^ —c{x,y) together with (a, vr) ^ (— a, — vr) . (2.23) 

(Note that under these reflections we also have a{x, y) —d{x, y), as we can see from 
(|2.7|) .) The reflection (cr, vr) (— cr, — tt) just shifts both asymptotic chiral angles 6± 
by the same amount vr, and clearly does not change the physics. Since this reflection 
interchanges b{x, y) and c(x, y) without affecting the physics, these two objects must 
have the same singularities as functions of a;, consistent with isospectrality of Li, and 



^This is true for short range decaying potentials on the whole real line. For periodic potentials 
both operators may have that ui'^ — tt^ mode in their spectrum Strictly speaking, (to the 

best of our knowledge) only the case it — Q appears in the literature on supersymmetric quantum 
mechanics. 



13 



2.2 The Diagonal Resolvent 

Following ||2y, |21| we define the diagonal resolvent {x\iD^^\x) symmetrically as 

/ A{x) B{x) 

^ C{x) D{x) 



[x 



iD-^\x) 



- lim 

2 e^o+ 



/ a{x,y) + a{y,x) b{x,y) + b{y, x) \ 
V c{x, y) + c{y, x) d{x, y) + d{y, x) ) 



(2.24) 



y=x+€. 

Here A{x) through D{x) stand for the entries of the diagonal resolvent, which follow- 
ing ( p.7|) and (|2.1(]|) have the compact representation^ 



B{x) 
C{x) 



bi{x)b2{x) 

Ci{x)c2{x) 



D{x) 
A{x) 



1 [d^ + 2a{x)] B {x) 

2 uj + 7r(x) 

i_ [d^ - 2(7{x)] C (x) 
2" 



(2.25) 



uj — 7r(x) 

We now use the generalized "supersymmetry" of the Dirac operator, which we 
discussed in the previous subsection, to deduce some important properties of the 
functions A{x) through D{x). 

From (|2.25| ) and from ( p.5|) we we have 

« 9^ - 2(T(a;) / CiC2\ i C2Q^ Ci + ciQ^ C2 



^^^^ 2 uj-n{x) [ Wj 2Wc cj-7r(x) 
Using ( |2.16| ) first, and then ( ^.151 ), we rewrite this expression as 



A{x) 



-(C261 + C162) 



i biQb2 + b2Qbi 



2VFc ' 2Wc cj + 7r(x) 

Then, using the fact that = (Eq. (|]2|)) and (12^251 ), we rewrite the last 
expression as 

A{x) 



i dx + 2a ( 6162' 



2 a; + 7r(x) V Wb 



1 {dx + 2a)B{x) 

2 uj + 7r(x) 



^A, B, C and D are obviously functions of uj as well. For notational simplicity we suppress their 
explicit Lo dependence. 
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Thus, finally, 

A{x) = D{x) . (2.26) 

Supersymmetry renders the diagonal elements A and D equal. 

Due to ( p.25| ), A = D is also a first order differential equation relating B and C . 
We can also relate the off diagonal elements B and C to each other more directly. 
From ( ^.25[ ) and from ( |2.15| ) we find 



After some algebra, and using ( p.20| ), we can rewrite this as 



h' h' 



[uo + TiyC = a'B + aB' + 



The combination h'^h'2/Wi, appears in B" = (6162/^^6)"- After using Lhbi^2 = to 
eliminate 6'/ and 62 from B", we find 

= -5" - + 7T^ - a' -uj^ ^ B 

Thus, finally, we have 

- (a; + n^C = -B" +{a- , \ B' - ( - a' - uj^ + B . (2.28) 

In a similar manner we can prove that 

{uJ - nfB = -Ic" +(cr- C + (n' + a' - + C . (2.29) 

We can simphfy ( p^.28| ) and ( ^.29|) further. After some algebra, and using ( ^.25| ) we 
arrive at 

C{x) = —^d^Mx)-^^^^B{x) 

B{x) = ' d,A{x)-'^^^^C{x). (2.30) 

UO — -K^X) UO — TX{X) 

Supersymmetry, namely, isospectrality of and Lc, enables us to relate the diagonal 
resolvents of these operators, B and C, to each other. 
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Thus, we can use ( |2.25| ), ( p.26| ) and ( |2.30| ) to eliminate three of the entries of the 
diagonal resolvent in ( p.25|) , in terms of the fourth. 

Note that the two relations in ( ^.301 ) transform into each other under 



B ^ —C simultaneously with (a, vr) 



-a, -7r) 



(2.31) 



in consistency with (|2.23|) . The relations in ( |2.30D are linear and homogeneous, with 
coefficients that for dxTT{x) 7^ do not introduce additional singularities in the u 
plane. Thus, we see, once more, that B and C have the same singularities in the u 
plane. We refer the reader to Section 4 in for concrete examples of such resolvents. 

The case tt{x) = brings us back to the GN model. In the GN model, our B 
and C, coincide, respectively, with uR- and —ljR+, defined in Eqs. (9) and (10) in 
[|19|. With these identifications, the relation A = D (Eq. ( p.26| )) coincides essentially 
with Eq. (18) of p|. The relations (^^) and were not discussed in but 

one can verify them, for example, for the resolvents corresponding to the kink case 
a{x) = mtanhmx (Eq. (29) in [l^), for which 

c = — , : B= I I -1 c. 



2\Jvn? — up' ' 



UJ 
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2.3 Bilinear Fermion Condensates and Vanishing of the Spa- 
tial Fermion Current 

Following basic principles of quantum field theory, we may write the most generic 
flavor- singlet bilinear fermion condensate in our static background as 



{^aait,x)Taf3^Paf3it,x)),eg = N / ^ 



r(x| 



—I 



\x) 



reg 



A{x) B{x) 
C{x) D{x) 




(2.32) 



where we have used ( p.24|) . Here a = 1, ■ ■ ■ , is a flavor index, and the trace is taken 
over Dirac indices a,j3. As usual, we regularized this condensate by subtracting from 
it a short distance divergent piece embodied here by the diagonal resolvent 

A B \ ^ / imcosO uj + msin6' 

'vac 



{x \ — iD ^\x ), 



C D 



VAC 



-uj + msin6' imcosO 



(2.33) 



of the Dirac operator in a vacuum configuration <7y^^ = mcos9 and vr^^^ = msin6'. 



In our convention for 7 matrices ( |2.3|) we have 
/ A{x) B{x] 



A{x)+D{x) _^^ A{x)-D{x) 1 ^ ^. g(x) - C{x) p ^ B{x) + C{x) 
\ C{x) D{x) ' 2 2« 2 2 

(2.34) 

An important condensate is the expectation value of the fermion current (j^(x)). 
In particular, consider its spatial component. In our static background ((7(x), 7r(x)), 
it must, of course, vanish identically 



{j\x)) = 0. 



(2.35) 



Thus, substituting F = 7^ in (|2.32|) and using ( p.34| ) we find 



{j\x))=zN I ^[A{x)-D{x)]. 



(2.36) 
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But we have already proved that A{x) = D{x) in any static background {a{x),iT{x)) 
(Eq. (|2.26|) ). Thus, each frequency component of (j^) vanishes separately, and (|2.35|) 
holds identically. It is remarkable that the generalized supersymmetry of the Dirac 
operator guarantees the consistency of any static {a{x),TT{x)) background. 

We discussed < j^{x) >= in |]21[]. However, that analysis was incomplete as it 
considered only the asymptotic behavior of < j^{x) >, which misled us to draw an 
overrestrictive necessary consistency condition on the background. 

Expressions for other bilinear condensates may be derived in a similar manner 
to the derivation of < j^{x) > (here we write the unsubtracted quantities). Thus, 
substituting T = 7° in ( p32|) and using ( p^lj) , and ( p30| ), we find that the 

fermion density is 

{j\x)) = ^N — [B{x) - C{x)] = ^N — ^ ' • (2.37) 

J Ztt j Ztx uo -\- T{\X) 

Similarly, the scalar and pseudoscalar condensates are 

(^(a;)^(x)) = iV / ^ [A{x) + D{x)] = 2iV | ^ D{x) , (2.38) 

and 

(^(a;)7'^(x)) =N — [B{x) + C{x)] =N — ^ ' ■ (2-39) 

J Ztt J Ztt lo + 7T[x) 
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3 The Saddle point Equations and Reflectionless 
Backgrounds 



For static backgrounds (o"(x), 7r(a;)) we have the (divergent) formal relation 
1 I , r doj . 1 , . 



duj 
2^ 



A{x) B{x) 
C{x) D{x) 

(see Eq. ( |2.32| )). Therefore, using ( |2.3| ) and ( p.26[ ), the bare saddle point equations 
(|1.5|) for static bags are 

HWJy . = -^-iv/^IBW + CMl^o, (3.1) 

I stjOitic 

where (cr(x), t^{x)) are subjected to the asymptotic boundary conditions (2^). As was 
already mentioned following ( |2.2| ), we further assume that <j{x) and 7r(x) tend to their 
asymptotic boundary values on the vacuum manifold at an exponential rate which 
is determined, essentially, by the mass gap m of the model, as typical of sohtonic 
configurations. 

The u;- integrals in ( p.l|) are divergent. For bounded bag profiles which satisfy the 
boundary conditions (|2.2|), the diagonal resolvent ( |2.24| ) tends, for large to that 



of the vacuum background (|2.33|) . Thus, we note from (p.33|) , that while each of the 



integrals / duo B{x) and / dujC{x) diverges linearly with the ultraviolet cutoff, their 
sum diverges only logarithmically, as does / du A{x). The saddle point equations for 
vacuum condensates, i.e., the gap equations 



,2 / 2vr 



^vAc rduj 



9' 



J ^ [By.c + C.Ac] = , (3.2) 
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exhibit the same logarithmic divergence, of course. Thus, we can take care of the 
UV divergence in (|3.1|) by subtracting from these equations the corresponding gap 
equations. 

We now concentrate on the subtracted saddle point equation for cr(x) 



cr X —a 



VAC 



duo 



[A{x) - A 



VACl 



(3.3) 



2Ng'^ Jc 2tt 

The integration contour C in ( |3.3|) is commonly^ taken as indicated in Fig.(l), which 
shows qualitatively the spectrum of the Dirac equation in a bag background. 



Im(x) k 



m 



CO 



ReOd 



m 



Fig. 1: The contour C in the complex uj plane in Eq. (3.3). The continuum states appear as 
the two cuts along the real axis with branch points at ±to (wiggly lines), and the bound state is the 
pole at uJi. 



Besides the continuum states in that spectrum (the two cuts corresponding to the 
Fermi sea of negative energy states u < —m, and scattering states with u > m), 
there are bound states within the gap —m < uj < m, which trap fermions into the 
bag. One such bound state is indicated in Fig. (1) as the pole at u = Ui. The 
detailed calculation to determine bound state energies like Ui is discussed in p|, ^ 
(after establishing the reflectionless property of the background). We stress that cr^^c 
in ( p.3|) can be the a component of any point on the vacuum manifold cr^ + vr^ = 



(and similarly for vr^^^, which appears in the subtracted equation for 7r(x)). 

^See e.g., Section 4 of |Q for a detailed discussion. If the bound state at energy wi traps ny 
fermions, then routing the contour C as indicated in Fig. (1) means that we are discussing the sector 
of the model with fermion number Uf. 
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The contour integral in ( |3.3| ) is most conveniently calculated by deforming the 
contour C into the contour C shown in Fig. (2). (This is allowed, since the subtracted 
integral in ( |3.3| ) converges.) The "hairpin" wing of C picks up the contribution of 
the filled Fermi sea, and the little circle around the simple pole at u = ui is the 
contribution of fermions populating the bound state of the "bag" . 



/m 0) A 



- m 



CO 







m 



Fig. 2: The deformed integration contour C' used in Eq. (3.3). 



^^=^ ReOd 



Let us study the spatial asymptotic behavior of ( |3.3| ). From ( |A.13| ) we have 

ik — a(—oo) 



2k 



ik + cr(+oo) 



■r2{k)e 



-2ikx 



ri{k)e 



2ikx 



ri{k)e'"'^ + 0{e 



-const. 'mx\ 



(3.4) 



where a(±oo) are the appropriate vacuum boundary values of a from ( |2.2|) , and 
f'i{k),r2{k) are the reflection coefficients defined in (|A.1|) . 

Thus, for example, studying ( p.3| ) as s ^ oo, we see from ( p. 4] ) that 

cr(x) — cr(+oo) r du ik + (t(+oo) 

x^+oo JFermi sea 271 



(3.5) 



2Ng'^ x^+oo JFermi sea 271 2k 

The first term in ( |3.5| ) is the contribution coming from the Fermi sea (i.e., the "hair- 
pin" wing of C). The second, exponentially small term on the right hand side of 
( |3.5| ) comes from the bound state pole (i.e., it is proportional to the bound state wave 
function squared). Due to the asymptotic boundary conditions on a{x), the left hand 
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side of ( |3.5| ) is also exponentially small as a; oo. Thus, the first term on the right 
hand side of ( |3.5| ) must have an exponentially small bound as x — > oo. 



We now change the variable to k = \J lJ^ — m^. When mapping into the fc-plane, 
the lower wing of the cut in Fig. (2) is transformed into k = l/^le**^, and the upper 
wing is transformed into —\k\e~''^, with 6 0+Q. Thus, we may write the dispersion 
integral on the right hand side of (|3.5| ) coming from states in the Fermi sea as — A(x), 
where 



A(x) 



dk {ik + (t(oo)) ri{k)e^'^'' + {k 

4:71 



-k) 



Re 



dk {ik + a {oo)) ri{k)e^'''^'^ 
2^ ~ 



(3.6) 



where in the last equality we used the reflection property ri(— A;) = rl{k) (|A.2[) of 
ri{k). 

The function A{x) has to die off at least at an exponential rate as a; ^ oo. Thus, 
we are to study the asymptotic behavior of 



G{x) 



dk {ik + a {oo)) ri{k)e^'-^'^ 
2^ ~ 



(3.7) 



at large x. To this end we have to invoke some of the general properties of the 
reflection coefficient ri{k) of the operator L;, in ( |2.9| ). 



Due to the boundary conditions ( pl2|) on the background fields a{x) and vr(x), the 
operator Lf, tends exponentially fast (in x) to its asymptotic free particle form. Thus, 
its "scattering potential" is localized in a finite region in space. 

From the literature on scattering theory (in one space dimension) we know 
that the refiection coefficient r{k) of Schrodinger operators with short range potential 
wells[^ is analytic on the real k axis (and generally follows the threshold behavior 



"'^'^For example, just above the cut w 
(5^0+. Thus, just above the cut, k 



m 



m — \lo + m|e**^^ ^\ with 



^^We tacitly assume that the potential wells in question tend to the same asymptotic value at 
X = ±(X) (as Lb does with (t(x), 7r(2:) satisfying (^)), and that they do not have any barriers above 



these asymptotic values. 
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r{k) = — 1 + a/c + ■ ■ ■) and dies off like as oo, i.e., at large kinetic energy. 
Strictly speaking, the discussion of these issues in the various references in con- 
centrates mostly on Schrodinger operators of the standard form —d^ + V{x), but the 
arguments leading to the conclusions about the behavior of r{k) may be easily gen- 
eralized to the scattering theory of the Sturm-Liouville operators and in ( |2.9| ) 
(with cr(x) and tt{x) relaxing fast to ( p.2|) ). 

Therefore, in deriving the asymptotic behavior of G{x) and A(x) we may use as 
an input that ri{k) is analytic on the real k axis and that it decays at least as fast 
as 1/k as k oo. Given these properties of ri(fc), we are allowed to expand G{x) 
in powers of 1/x in the most natural way, namely, by repeatedly integrating by parts 
over k in (|3.7| ). 

Thus, for example, after three integrations we find 

(T(oo)ri(O) imr i{0) + ma {oo)r[{0) 

iTTLriX) = 1 —: 

zimx [2imx)'^ 

-(T(oo)ri(O) + 2im'^r[{{]) + m'^a{ooyi{{)) 
{2imxY 

dk- -— ^ ' riik) , 3.8 

and so on and so forth. Clearly, the remaining integral in each step is subdominant by 
a power of 1/x relative to its predecessor, and thus, the expansion of G{x) generated 
in this way is an asymptotic expansion. 

From (p.8|) (or by working out a few more terms in the asymptotic expansion if 
necessary) the following pattern emerges: the coefficient of (l/2ix)^""^^ is a linear 
combination of the form 

2n 

i=o 

with real coefficients cj, and the coefficient of (l/2zx)^" is a linear combination of the 
form 



2n-l 

i 

j=0 
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with some other real coefficients cj. 

From the reflection property ( |A.2| ) ri{—k) = rl{k) we immediately conclude that 
r["^(0) = {—l)^r["\o)*: the even derivatives r^''\o) are real, and the odd derivatives 
^(2j+i) ^g-j imaginary. Thus, the linear combinations J2 CjVvi^ (0) are real, which 
makes each term on the right hand side of the asymptotic expansion ( |3.8| ) of G{x) 
pure imaginary. 

Using this result in ( |3.6|) we conclude that all terms in the asymptotic expansion 
of A(x) in powers of 1/x vanish. Thus, A(x) vanishes faster than any power of 1/x 
as X ^ oo. This is consistent with our expectation that A (a;) vanishes at least at an 
exponential rate when x oo. 

This concludes our discussion of the subtracted saddle point equation ( |3.3|) for 
(t(x) and its asymptotic behavior. 

We can repeat the same story for the subtracted saddle point equation for •n^x) 

^i^J""" = - / ^ - B,,a) + {C{x) - C,,,)] . (3.9) 

In a similar manner to our derivation of ( p.5| ) and ( |3.6| ), we can show that 

7r(x) — 7r(oo) . „ „ .1 . r 1 1 
-—: — > exponentially small contribution of bound states 

+ — Trfcxo)) — (Vk'^ + + n(oo)) '^[°°}~^^.'! 

(3.10) 



-Re 



dk{V^T7^-7r{cx^))-{VW+7^ + n{oo)r^0±^^ ^^^^ 

-^=^= — — ri{k)e 

271 V A;^ + 



Due to the asymptotic boundary conditions on tt{x), the left hand side of ( |3.10[ ) is 
exponentially small as x oo, which bounds the dispersion integral on the right 
hand side of (|3.10|) . As in our analysis of ( p. 71) , we expand the integral in the square 
brackets in ( p.lO|) in powers of and similarly to (|3.8| ), we can show that all terms 



in that asymptotic series are pure imaginary. Thus, the right hand side of (|3.1CI|) 
vanishes faster than any power of 1/x, consistent with the boundary conditions on 
7r(x). 

We conclude that the asymptotic behavior of the static saddle point equations 
(|3.3|) and ( p.9|) is consistent with the asymptotic boundary conditions ( |2.2| ) on a{x) 
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and tt{x) for any reflection amplitude ri(fc): all terms in the asymptotic expansions 
of the dispersion integral in the square brackets in ( p.lO| ) and also of G{x) in ( p.7| ) in 
powers of 1/x are imaginary. 

Contrary to the argument we made in pT| , the reflectionless property of the solu- 
tions (j{x) and tt{x) of does not emerge as a necessary condition from consistency 
of the asymptotic behavior of (|3.5| ) and (|3.1CI|) and the boundary conditions on the 
background fileds. 
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Appendix: Asymptotics of the Dirac Resolvent 



In this Appendix we discuss the spatial asymptotic behavior of the diagonal resolvent 
of the Dirac operator ( |2.24| ). 



According to our discussion in Section 2.2, given, for example, B{x), we may 
determine -D(x), A{x) and C{x) using ( p.25| ) first, then ( |2.26| ) and finally, ( p.28| ). 

Thus, it is enough to determine the asymptotic behavior of B{x). According to 
(|2.25|) , B{x) = bi{x)b2{x)/Wb, so we need to determine the asymptotic behavior of 
61 (x), 62(3;), the fundamental solutions of Lbb{x) = 0. 

By definition, according to (|2.11|) , and with the standard choice ( p.l7 ) of the 



coefficients A^^ = A^^^ = 1, these functions satisfy 



bi{x) — > e '^'^ , b2{x) — ^ e*''^ 
x— >— 00 



at the two opposite sides of the one dimensional world. Thus, what remains is to 
determine the asymptotic behavior of each of these functions on the other side of the 
world. Since a and vr relax to the vacuum manifold as x — ±cxd (Eq. (p.2|) ), the 
operators Lf, and Lc in ( p?9|) degenerate into free particle operators, and thus we must 
have 

^ ' ti{k) ti{k) 

b^ix) -^rre""" + '^e-"'\ (A.l) 

where ti, 2(^)5 ''"i,2(fc) are the appropriate transmission and reflection amplitudes, re- 
spectively. 

It is enough to consider /c > 0, since the operators and L^. in ( |2.9D are real and 
thus bi^2{x, —k) = bl 2(2;, k), leading to 

ti,2{-k) = tl^{k) and ri,2(-fc) =r*2(fc). (A.2) 

Thus, bi{x) corresponds to a setting with a source at x = +00 which emits to the 
left, and 62(2;) describes a source at x = ~oo which emits to the right. 
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The wronskian of 61 (x) and 62 (a^) (Eq. ( |2.13| )) 

b2{x)b[{x) - bi{x)b'2{x) 

Wb{uj) = ■ — 

u + n{x) 

is independent of x. Thus, evaluating it at x ^ ±00 we find 

^ t2{k){uj + n{-oo)) ^ ti(A;)(cu + 7r(+oo)) ' ^ ' ' 

and thus, 

ti{k){uj + 7r(+oo)) = t2ik){uj + 7r(-oo)) . (A.4) 

Like the wronskian of {bi{x),b2{x)}, the wronkians of the pairs of independent 
solutions {bi{x) , bl{x)} and {62(3^), &2(^)} (here we assume that k is real) are also 
independent of x. In fact, these wronskians are proportional to the Schrodinger 
probability currents carried by bi{x) and 62 (3^)5 respectively. Thus, using ( |2.11|) and 
(|A. 1|) to evaluate each of these wronskians at a; = ±00 and then equating the results, 
we obtain 

bl{x)b'i{x) - (x) 



W,[bubl] 



uj + 7r(x) 
-2ik -2ik l-lriP 
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Wh[b2,b;] 



u; + 7r(— 00) + 7r(+oo) \t 

b*2{x)b'2{x) - b2{x)b*^ (x) 
UJ + 7r(x) 

2ik 2ik 1 - |r2p 



UJ + 7r(+cxD) UJ + 7r(— cxd) \t 
Thus, we deduce the generahzed unitarity conditions 



12 



(A.5) 



1 - 




ri 


2 




2 


UJ + 7r(+oo) 




1^1 


2 


1 - 




r2\ 


p ti; + 7r(— 00) 



(A.6) 



The usual unitarity condition + |rp = 1 holds for backgrounds in which 7r(+oo) 
7r(— 00). 
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Putting every thing together, we finally learn that 

1 + r2(fc)e-2*'=^ 



B(x) 



B(x) 



-2ik 



u + 7r(— oo)) 



as-— ++CXD 



1 + ri{k)e 
-2ik 



2ikx 



-{uj + 7r(+oo)) 



(A.7) 



Recall from ( p.33|) that 



UJ + 7r(±oo) io + 7r(±oo) 



-2ik 



(A.8) 



corresponds to the appropriate vacuum configurations ((7(±oo), 7r(±oo)). Thus, we 
can rewrite ( |A.7| ) as 



B(x) 



-2ikx 



B{x) 5,,,(+oo)(l + ri(fc)e2^^-) 

z— >+oo ^ ' 



(A.9) 



Using ( |2.26| ), ( p.27[ ) and ( |2.28| ) we then find the asymptotic behaviors of the re- 
maining entries of the diagonal resolvent ( p.24| ): 

ik — a(—oo) 



and 



A{x) = D{x) 



A{x) = D{x) 



DvAci-^) + 



2k 



■r2{k)e 



-2ikx 



X—f + OO 



I . ik + cr(+oo) , 2,-. 



C{x) 



C{x) 



x— >+oo 



cr( — oo) — ik _2ii,x 
a(-oo) +ik ^ ' 



(t(+cx)) +ik , 2,-t.r 
1 + — ^- '- -ri{k)e^"''' 



with 



from (|;33D. 



(j(+oo) — ik 
ijj + 7r(±oo) ijj — 7r(±oo) 



2^/rn? — uj"^ 



2ik 



(A.IO) 



(A.ll) 



(A. 12) 
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Finally, we can write these results more compactly as 



/ A{x) B{x) 
V C{x) D{x) 



a;— >+oo 



/ A{+oo) B{+oo) 
V C{+oo) D{+oo) 



+ 



VAC 



I i/c+(T(+oo) 



\ ^ V AcV^'-^ ) a{:^oo)-%k -2k I 



ri{k)e 



2ikx 



A{x) B{x) 
C{x) D{x) 



X—f — OO 



/ A{-oo) B{-oo) 
V C(-oo) D{-oo) 



+ 



VAC 



I ik—a{~oo) 
2k 



\ ^VAci ^)a{-oo)-ik 



cr(—oo)+ik ik—cr(—oo) 
2k 



r2{k)e 



-2ikx 



(A.13) 
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